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ABSTRACT 

A special purpose Finite element composite analysis program (FECAP) for analyzing composite 
material behavior with a microcomputer is described. The formulation assumes a state of 
generalized plane strain in a material consisting of two or more orthotropic phases. Loading 
can be mechanical and/or thermal. The theoretical background, computer implementation, and 
program users guide are described in detail. A sample problem is solved showing the required 
user input and computer generated output. 


NOMENCLATURE 


Matrix Notation: 

| | rectangular or square matrix 

| matrix inverse 

[ ]T matrix transpose 

{ ) column vector 

{ } T row vector (transpose of column vector) 
Symbol Definitions: 

All symbols are defined in the text. 


INTRODUCTION 

Advanced composite materials have gained wide use in the aerospace industry over the last 20 
years because of their high specific strength and stiffness, and low coefficient of thermal 
expansion (CTE). Design of composite structures requires the analysis of composite material 
behavior. Composite materials, for the purpose of this report, are defined as the combination 
of at least two distinct materials to form one nonhomogeneous anisotropic material. The 
anisotropic nonhomogeneous nature of composites has limited the number of closed-form 
analytical solutions available to analyze their behavior. Numerical techniques have been widely 
applied to many problems in composite material behavior. The finite element method has been 
one of the most popular of the numerical techniques. 

Finite element analyses generally require large amounts of computer memory, and therefore in 
the past have been mainly used on large main-frame computer systems. However, recent 
advances in microcomputer technology have made finite element analysis on desk-top machines 
practical. Many large general purpose finite element programs are now commercially available 
for microcomputers. This report will describe a special purpose finite element composite 
analysis program (FECAP) developed at NASA-LaRC to analyze composite material behavior 
on a microcomputer. The theoretical formulation is presented first, followed by a description 
of the computer program and a users guide. A sample problem is solved showing the required 
user input and computer generated output. 
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Identification of commercial products in this report is provided to adequately describe the 
products and does not constitute official endorsement, expressed or implied, of such products 
or manufacturers by NASA. 


THEORETICAL BACKGROUND 


Problem Description 

The problem under consideration is an orthotropic body, consisting of one or more orthotropic 
materials, in a state of generalized plane strain. The material is assumed to be linear elastic and 
temperature independent. The coordinate system, and the different levels of analysis possible 
for a continuous fiber reinforced composite (i.e. laminate or constituent level), for this problem 
are shown in fig. 1. The displacement field for this problem may be written as: 


u(x,y,z) = U(y,z) + (e xo )x 

v(x,y,z) = V(y,z) (1) 

w(x,y,z) = W(y,z) 

where u, v, and w are the total displacements in the x,y, and z directions, respectively, and U, 
V, and W are the unknown displacement functions to be determined. e x0 is a constant 
uniform strain in the x direction, and may be treated as either a known or unknown quantity. 
The displacement functions U, V, and W are only functions of y and z. Therefore the problem 
becomes two dimensional even though all three components of displacement, and all six stress 
and strain components are present. This type of problem is often referred to as quasi-three 
dimensional. 


Finite Element Geometry 


The finite element formulation [1] is presented for both 4-node and 8-node isoparametric 
elements. The global and local coordinate systems for this type of element are shown in fig. 2. 
The relationship between the global and local coordinates is written as 


y = Njyi i=l,4 or 1,8 
z = Njzi i=l,4orl,8 


(2) 


where the interpolation functions Nj are functions of this local coordinates (i;,ri) , and yj and zj 
are the global coordinates of the nodes. 

The interpolation functions for the 4-node element are linear, and are defined as 

Ni = (l/4)(1+^ 0 )(1+Ti 0 ),i=l,4 (3) 


where = (^)(%i)» "Ho = 01)010. and ^ and T| range from -1 to +1 at the comers. The 
interpolation functions for the 8-node elements are quadratic, and are defined as 

Corner nodes - 

Ni = (l/4)(l+^o)(l+Tio)(^o+Tio-l), i=l,3,5, and 7 (4) 

Mid-side nodes - 

Nj = (1/2)(1-^)(1 +ti 0 ), 5i=0, i=2 and 6 (5) 
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Nj = (l/2)(l+5o)(l-‘n 2 ), UpO, i=4 and 8 

Both and r| 0 have the same definition as given for the 4-node interpolation functions. 


( 6 ) 


Isoparametric elements use the same interpolation functions for the geometry and the nodal 
variables (i.e. displacements). Therefore, the unknown displacement functions of equ. (1) 
may be written as 


U(y,z) = N iUi 

V(y,z) = Njvj (7) 

W(y,z) = NjWi 

where uj, vj, and wj are the unknown displacements at the nodes. 

Finite Element Formulation 

The finite element formulation used in FECAP is displacement based and requires the 
minimization of the total potential energy, PE, for each element with respect to the unknown 
variables. This procedure leads to a set of linear simultaneous equations relating the unknown 
nodal displacements to the applied loads. The equations for each element are assembled into a 
global system, the boundary conditions are applied, and the system is solved for the nodal 
displacements. A detailed description of the entire finite element formulation is given in 
reference [1]. A brief description is given in this report. 

The total PE of an element is the sum of the strain energy, Us, and the work of the external 
loads, Wl. The strain energy for an element is written as 

U S = l/2Jvol (<3> {e m } T dvol (8) 

where (o) and {e m } are the element stresses and mechanical strains, respectively, and have 
the form 


{<j)^ - {o x , Oy, o z , ty Z , T Z X’ ^xy}'^' (9) 

( e m } T =( e mx> e my > E mz>Ymyz»Ymzx>Ymxy ) ^ (10) 

The work of the external loads is written as 

W L = -{q}{f} T -{ex 0 }(F X0 ) (11) 

where (q) are the nodal displacements, { f) are the applied nodal loads, and F xo is the average 
axial force. 

The total strain for a linear elastic system is written as the sum of the mechanical strain and the 
free thermal strain, or 


(el = {e m } + ( e t) 

Rearranging to solve for the mechanical strain leads to 


( 12 ) 
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{e m ) = { £ } ' Ut) 


(13) 


Minimization of the total PE with respect to the nodal displacements requires that (o) and 
(e m ) in equ. (8) be expressed in terms of the displacements, and is accomplished by using the 
strain-displacement relationships of linear strain theory. These relationships are written in 
matrix form as 


{£} = [L]{u} (14) 

where { u } is the vector of total displacements, and the matrix [L] is defined as 


d/dx 

0 

o - 

0 

d/dy 

0 

0 

0 

d/dz 

0 

d/dz 

d/dy 

d/dz 

0 

d/dx 

d/dy 

d/dx 

0 - 


Note that engineering shear strain has been used in the above relationships. 

Combining equs. (1) and (7) and substituting into equ. (14) leads to an expression for strain in 
terms of the nodal displacements written as 

{e}=[B]{q} + {e xo ) (16) 

where 

{q} T = {ui.vi.wj, ... u n ,v n ,w n ) (17) 

and 

{£xo) T = {e xo .0,0, 0,0,0) (18) 

The matrix [BJ (6 x 3n) is defined as the matrix product of [L] and a partitioned matrix 
involving the interpolation functions written as: 

[B] = [L] [LI](N ! )» [I](N 2 )| ... [I](N n )] (19) 

where I is a 3x3 identity matrix. The subscript n in equs. (17-19) takes on values of 4 or 8, for 
the 4-node or 8-node elements, respectively. Substitution of equ. (16) into equ. (13) leads to 
an expression for the mechanical strains in terms of the displacements, written as 

( £ m } = [B] {q} + {e xo } - (e t ) (20) 

where the free thermal strain vector is defined as 

{e t ) = {a}(AT) (21) 
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and AT is a uniform change in temperature. 


The stresses are expressed in terms of the displacements through the material constitutive 
equations, written in general form as 


{a} = [Cbar]{e m } (23) 

The specific form of [Cbar] will be discussed in the next section. 

Equations (20) and (23) are substituted into equ. (8) to write an expression for the strain 
energy in terms of the displacements as 

Us = ltfJvol (([B]{q)+{ex 0 }-{£ t ))T[Cbar] 

([B]{q}+{e xo }-Ut})) dvo1 (24) 

Two cases must be considered for the minimization of the total PE (equs. (1 1) and (24)) with 
respect to the unknowns. First, for the case of a known e x0 , the PE is minimized with respect 
to only the nodal displacements, (q). This minimization results in an expression written as 

[K]{qJ = {F} (25) 

where 

[K] = J vo i [B]T[Cbar][B]dvol (26) 

and 

(F) = -J vol [B]T[Cbar]({e xo }-{e t })dvol + {f} (27) 


The case of an unknown e x0 requires that the PE be minimized with respect to both the nodal 
displacements, {q}, and e x0 . This results in an expression of the general form 


[K*]{q*J = {F*} 

which may be expanded and written in partitioned form as 


[K] ;{K x n 

i 

(q) 


(F) 

(k x } t ; k xo 

e xo 

" 

X 

o 


where [K], {q } , and {F} have the same definitions as previously given, and 


(28) 


(29) 


{K x } = l’st column of 


lvollB] T | Cbarjdvol 


(30) 



and 


Kxo - 


(31) 


F x = l'st element of J vo i[Cbar](e t }dvol (32) 

Two steps are required to compute the elements of either [K] or [K*] and {F} or {F*}. First, 
the derivatives of the shape functions, Nj, with respect to the global coordinates, (y,z), must be 
evaluated for terms involving [B]. However, the shape functions are in terms of the local 
coordinates, (£,t]). The derivative evaluation is acomplished by use of the Jacobian matrix in 
an expression written as 


3Nj/3y 


aNi/a^ 


= [JJ- 1 


5Nj/dz 


dNi/dr| 


where the Jacobian matrix is defined as 


“aNi/as ...dN n /aC 


yi z i 

_3N 1 /Dti ...3N n /0T}_ 


_yn z n- 


(33) 


(34) 


The subscript n in equs. (33) and (34) takes on values of 4 and 8 for the 4-node and 8-node 
elements, respectively. 

The second step is to numerically evaluate the integrals in equs. (26-32). This is accomplished 
by first transforming the integral into the (^,T|) coordinate system by noting that 


dvol = |J I (d^(dn) 


(35) 


where Lf I is the determinate of the Jacobian matrix. The limits of integration become -1 to +1 

on both £, and rj. Gauss quadrature is used for the integration and all of the integrals can be 
written in the form 


n n 

Integral = Z X HjHjG(£i,r|j) Ll I (36) 

i=l j=l 

where Hj and Hj are the weight functions and G(^,ri) is the function to be integrated. A 2x2 
Gauss quadrature (n=2) is used for both the 4-node and 8-node elements. The weight factors 
and coordinates of the Gauss points are given in [1]. 

Equations (25) and (29) represent a system of linear simultaneous equations for one element. 
These element equations are assembled into a global system of equations by requiring 
continuity of the displacements at the nodes of adjoining elements. The prescribed boundary 
conditions are imposed and the resulting system of simultaneous linear equations is solved for 
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the unknown nodal dispalcements, and if necessary e xo . The element strains and stresses are 
computed from the nodal displacements by using equs. (14) and (23). The strain energy 
density (SED) is computed from 

SED =[a x (£ x -a x AT)+Gy(£y-ayAT)+a z (e z -a z AT)+Gy Z (Yy Z - a yz AT) (37) 


+<7 zx(Yzx' a zxAT)+G X y(Y X y-a X yAT)]/2 

Resultant internal nodal forces [1] are computed from 

(P) = Jvol [B]T{cr)dvol (38) 

where { P) is the vector of resultant nodal forces for an element. The integration in equ. (38) is 
evaluated numerically by rewriting equ. (38) in the form of equ. (36). The summation of these 
nodal forces in the x, y, and z directions, respectively, must equal zero for static equilibrium, 
and is used as a check on the validity of the solution. 

It should be noted that the formulation presented above (i.e. minimization of the total potential 
energy) insures that IK) and [K*] will be symmetric and positive definite, after the imposition 
of boundary conditions which prohibit rigid-body motion. Therefore the inverse of [K] and/or 
[K*] exists and a solution can be found. 


Material Constitutive Equations 

A general form of the material constitutive relations (equ. (23)) was given in the preceding 
section. The specific form of these relations depends upon the relative orientation of the 
principal material axes with respect to the global axes. FECAP can analyze orthotropic 
materials (i.e. materials with 3 principal planes of symmetry) with two types of orientations 
relative to the global axe's. These orientations are refered to as plane and cylindrical, and are 
shown in fig. 3. The plane orientation involves a rotation about the z or 3 axis, and the 
cylindrical orientation involves a rotation about the x or 1 axis. The elements of [Cbar] for 
each orientation are given in the Appendix A. 

The angle of rotation, 0, for the cylindrical orientation varies with position in the y-z plane, and 
therefore varies within a given element. FECAP accounts for this variation by computing theta 
from the global coordinates at each of the Gauss points used for the numerical integration in 
equs. (26-27) and equs. (30-32). FECAP can also analyze materials having a combination of 
both orientations. 

The matrix {a} also depends upon the relative orientation of the principal material axes with 

respect to the global axes. The specific forms of {a} for the two orientations described above 
are also given in Appendix A. 


Solution Procedure 

As discussed above, the solution of equs. (25) or (29) first requires the imposition of boundary 
conditions on the system. The three types of boundary conditions allowed in FECAP are 
prescribed nodal displacements, applied nodal forces, and constrained nodal displacements. 

The first two involve standard procedures for appropriately modifying the system stiffness 
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matrix ([K] or [K ]) and load vector ({F} or {F }), and are well documented in texts on finite 
element analysis [1]. The application of constrained nodal displacements is not as well 
documented and requires further explanation. 

Constrained nodal displacements refer to a series of nodes whose displacements are 
proportional to one another but whose absolute values are unknown (e.g. v\ = V 3 = 4v£). 

This effectively reduces the total number of unknowns in the system. The system stiffness 
matrix and load vector are modified by combining the contributions of the affected nodes into a 
single row and/or column. These modifications do not destroy the symmetry of the stiffness 
matrix, but can increase the bandwidth (to be discussed subsequently). The specifics of the 
necessary modifications are given in reference [ 2 J. 

Determining the nodal displacements from equs. (25) or (29) requires finding the solution to a 
system of simultaneous linear algebraic equations. As stated earlier, [K] and [K*] will be 
symmetric and positive definite, and therefore an inverse will exist, and a solution is possible. 
The matrix [K] will also be banded (fig. 4) due to the assembly procedure in which a nodal 
displacement affects only those elements adjoining that node. Numerous solution algorithms 
exist for banded symmetric systems. These algorithms require storage of only the terms that lie 
on or above the diagonal and within the band, thus greatly reducing the amount of computer 
memory needed for the solution. However, [K*] is not banded due to the fact that every 
element contributes to e x0 . This type of matrix is sometimes refered to as an "arrowhead" 
matrix, and is shown schematically in fig. 4. A special solution algorithm was developed for 
FECAP that allows a banded storage scheme to also be used for [K*]. This is accomplished 
by storing the terms of [K*] which contribute to e xo as a separate vector. A standard triangular 
([L][Uj) decomposition ([ID is then used to solve the banded system for the unknowns, with 
terms from this separate vector used when needed. 


COMPUTER IMPLEMENTATION 


System Requirements 

FECAP uses the finite element formulation presented in the preceding section. FECAP is 
written in Hewlett Packard BASIC 3.0, and can run on any Hewlett Packard 9000 Series 200 
microcomputer with the BASIC 3.0 operating system. Although no special hardware or 
software enhancements are necessary, certain enhancements greatly increase program capability 
and performance, and are strongly recommended. The suggested system configuration is 
shown in Table 1. 

The suggested size RAM in Table 1. is 2 Mbytes, and controls the size of problems that can be 
solved with FECAP. An in-core solution scheme is used in the program and the size of the 
stiffness matrix, [K], is the major factor in determining the size of the problem. The amount 
of memory needed to store [K] can be computed from 

[K] Memory(bytes) = ((n)(3)(SBW))8 (39) 

where n is the number of nodes, 3 is the number of variables per node, SBW is the semi- 
bandwidth (fig. 4), and 8 is the number of bytes per word. For example, a medium size 
problem with 300 nodes and a semi-bandwidth of 75 requires 540 Kbytes. Loading the 
BASIC 3.0 operating system, and loading and compiling FECAP requires an additional 720 
Kbytes. Therefore, the total RAM required for execution of this example problem is 1.26 
Mbytes. The suggested RAM size of 2 Mbytes will provide enough memory for problems of 
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this size and larger. It should be noted that various factors can affect the semi-bandwidth, such 
as node numbering and constrained nodal displacement boundary conditions, and thus affect 
the memory requirements for a given problem. 


Pre- and Post-Processimi 

Pre-processing for a finite element program involves generating the global nodal coordinate 
data and the element connectivity data for the finite element mesh. FECAP is written such that 
this data is read from a file that the user has previously generated. This file can be produced 
from relatively simple user written programs, or from one of the many commercially available 
mesh generation programs. Data files from commercially available programs will have to be 
modified so that the data is in the proper format required for FECAP. The required format is 
given in Appendix B. A very simple mesh generation program was written for the example 
problems presented in this report. 

Post-processing in finite element analysis involves computing secondary quantities such as 
stresses and strains, and then displaying these quantities in tabular and/or graphic form. 

FECAP computes the nodal stresses and strains for each element from the nodal displacements, 
and displays these results in tabular form. FECAP also has the capability to store these values 
in a data file. The format of this data file is given in Appendix B. This data can then be used in 
a contour plotting program for graphical display of results. An interactive contour plotting 
program that is easily modified to work with FECAP can be found in reference [3]. 


User's Guide 


FECAP runs interactively, and thus requires responses from the user during execution. A 
description of the computer prompts and appropriate user replies are given below. This 
description assumes a working knowledge of the HP BASIC 3.0 operating system. 

Load file FECAP (compile if compiler is available) and begin execution by pressing the RUN 
key. FECAP will begin execution with the data input softkey menu appearing at the bottom of 
the screen and a computer prompt that reads: 

DATA INPUT - Select a key 

The user responds by pressing on of the seven defined softkeys. Softkeys K0 through K4 
must be selected at least once and in any order, but must be selected before K5 or K6. 

Softkeys K0 through K4 may be selected more than once if necessary to correct data. The 
selection of softkey K5 is optional, but must occur before the selection of K6. Softkey K6 
must be the last softkey selected. A description of each of the softkeys is given below. All 
user responses are entered by pressing the ENTER key. Multiple entries in the user response 
are seperated by commas. 

Title (K0): 

Computer prompt - Input Problem Title, (<= 80 char.) 

User response - Problem title, must be less than or equal to 80 
characters in length. 


Material Props (Kl): 

Computer prompt - Input Number of Materials 
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User response - Number of different materials, must be 
consistent with mesh data file. 

Computer prompt - Is a Cylindrical Material Orientation Required? (Y/N) 

User response - Is a cylindrical material orientation required for this problem, yes (Y) or 

no (N). Note, a cylindrical material orientation depends upon the mesh 
geometry and effects every element. Plane orientations can be used in 
conjunction with the cylindrical orientation. 

Computer prompt - For Material No. n (top of screen), Input Mat. Prop’s. 

E1,E2,E3,G23,G13,G12,NU23,NU13,NU12 
User response - Elastic orthotropic Young's moduli, shear moduli, and Poisson's ratios 

in principal material coordinates for n'th material. 

Computer prompt - For Material No. n (top of screen), Input Thermal Prop's. 

ALPHA 1,ALPHA2,ALPHA3 

User response - Orthotropic coefficients of thermal expansion in principal material 

coordinates for n'th material. 

Computer prompt - For Material No. n (top of screen). Input Number of Angles 

User response - Number of different plane orientation angles for n'th material, must be 

consistent with mesh data file. 

Computer prompt - For material No. n (top of screen), Input Angles in degrees (one at a 

time followed by ENTER key) 

User response - Plane orientation angles for n'th material. 

Note: The above four prompts are repeated n times, where n is the 
number of materials 

Mesh Data (K2): 

Computer prompt - Input Mesh Data File Name and MSUS (Ex.Fn, "700,0") 

User response - Mesh data file name and mass storage unit specifier (MSUS). The 

proper format for the mesh data file is given in Appendix B. 

Computer prompt - Input the Maximum Allowable Semi-Bandwidth (Must be <= Array 

dimensions) 

User response - Maximum allowable semi-bandwidth, must be consistent with array 

dimensions in main program. 

Boundary Data (K3): 

The boundary condition input softkey menu is displayed with the following prompt 
BOUNDARY CONDITION DATA - Select a Key 

These softkeys control how the boundary condition data is entered. Softkey K0 is for input 
from a previously defined file. Softkey K1 is for keyboard input. A description of each is 
given below. 

File (K0): 

Computer prompt - Input Boundary Condition File Name and MSUS (Ex. Fn, "700,0") 
User response - Boundary condition file name and mass storage unit specifier (MSUS). 

The proper format for the boundary condition data file is given in 
Appendix B. 
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Keyboard (Kl): 

The keyboard boundary condition softkey menu is displayed with the following prompt 
NODAL BOUNDARY CONDITIONS - Select a Key 

These softkeys control the keyboard input of the allowable types of nodal boundary conditions. 
Softkeys KO through K2 may be selected in any order, and more than once if corrections tire 
necessary. Softkey K3 is used to store the current boundary conditions in memory to a file. 
Softkey K4 returns the program to the data input softkey menu after completion of boundary 
condition input. A description of each softkey is given below. 

Specified Disp (KO): 

Computer prompt - Specified Nodal Displacements: Input Direction of Specified Nodal 

Displacement (1-X,2-Y,3-Z) (Input 0 for End of Specified Nodal 
Displacements) 

User response - Global direction of specified nodal displacement, 1 for x direction, 2 for 

y direction, and 3 for z direction. A response of 0 is used to end 
specified nodal displacement input. 

Computer prompt - Input Beg'ning node #, Ending node #, Increment, and Specified Value 
User response - Beginning node number, ending node number, node number increment, 

and specified nodal value for the direction specified in the previous 
prompt. 

Note: The above two prompts are repeated until a 0 is input for the direction at which point 
the program returns to the boundary condition softkey menu. 

Applied Force (Kl): 

Computer prompt - Applied Nodal Forces: Input Direction of Applied Nodal Force (I-X,2- 

Y,3-Z) (Input 0 for End of Applied Nodal Forces) 

User response - Global direction of applied nodal forces, 1 for x direction, 2 for y 

direction, and 3 for z direction. A response of 0 is used to end applied 
nodal force input. 

Computer prompt - Input Beg'ning node #, Ending node #, Increment, and Specified Value 
User response - Beginning node number, ending node number, node number increment, 

and specified nodal force for direction specified in previous prompt. 

Note: The above two prompts are repeated until a 0 is input for the direction at which point 
the program returns to the boundary condition softkey menu. 

Constraint Eqs (K2): 

Computer prompt - Constrained Nodal Displacements: Input Number of Different 

Constraint Conditions (Input 0 for No Constraint Conditions) 

User reponse - Number of different nodal constraint conditions. A response of 0 is 

used when there are no constraint conditions. 

Computer prompt - Input Direction of Constrained Nodal Displacement (1-X,2-Y,3-Z) for 

Constraint n 

User response - Direction of constrained nodal displacement for n'th constraint 

condition, 1 for x direction, 2 for y direction, and 3 for z direction. 

Computer prompt - Input Beg'ning node #, Ending node #, Increment, and Scale Factor 
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User response - Beginning node number, ending node number, node number increment, 

and constraint scale factor (i.e ratio of displacements) for the direction 
specified in the previous prompt. 


Computer prompt - More nodes for this direction ? (Y/N) 

User response - Are there additional constrained nodes for the direction specified above, 

yes (Y) or no (N). If the response is yes the previous prompt is 
repeated. 


Computer prompt - More nodes for this constraint condition ? (Y/N) 

User response - Are there additional node directions for the n'th constraint condition, yes 

(Y) or no (N). If the response is yes then execution returns to the 
constraint direction prompt. 


Note: The above four prompts are repeated n times, where n is the number of constraint 

conditions, after which the program returns to the boundary condition softkey menu. 


Save Data (K3): 

Computer prompt - Input Boundary Condition File Name and MSUS (Ex. Fn,"700,0") 
User response - Boundary condition file name and mass storage unit specifier (MSUS) 

for storing current boundary condition data. 


Exit (K4): 

This softkey returns the program to the data input softkey menu. The description of the data 
input softkey menu is continued below. 

Global Loads (K4): 

The global load softkey menu is displayed with the following prompt 
GLOBAL LOAD CONDITIONS - Select a Key 

Softkeys KO and K1 define the type of global loading for the problem, and only one of these 
softkeys should be selected. Softkey K2 returns the program to the data input menu. A 
description of each softkey is given below. 

Axial Strain (KO): 

Computer prompt - Input Uniform Axial Strain and Temperature Change 

User response - Uniform applied axial strain and uniform change in temperature. 

Axial Force (Kl): 

Computer prompt - Input Uniform Axial Force and Temperature Change 

User response - Uniform applied axial force and uniform change in temperature. 

Exit (K2): 

This softkey returns the program to the data input softkey menu, the description of which is 
continued below. 


Print Input (K5): 

Computer prompt - Input Printer Code for Output - CRT(l), PRINTER(701) 

User response - Printer code for output results, 1 for CRT display and 701 for printer at 

address 701. 

Computer prompt - Store Results on Output File ? (Y/N) 


User response - Store results in an output file for later use with post-processor program, 

yes (Y) or no (N). If response is no the following prompt is skipped. 

Computer prompt - Input File Name and MSUS for Output Results (Ex. Fn, "700,0") 

User response - File name and mass storage unit specifier (MSUS) for output results. 

Computer prompt - Print Nodal Coordinates ? (Y/N) 

User response - Print the global y and z coordinates of the nodes, yes (Y) or no (N). 
Computer prompt - Print Element Connectivity Data ? (Y/N) 

User response - Print the element connectivity data (i.e node numbers of each element), 

yes (Y) or no (N). 

Computer prompt - Print Boundary Condition Data ? (Y/N) 

User response - Print the nodal boundary condition data, yes (Y) or no (N). 

After the above prompt the following message appears on the screen 

Press Execute Softkey to Begin Execution 

Pressing the Execute softkey in the data input menu starts the finite element calculations. 
During execution the following message appears flashing on the screen 


***** program EXECUTING, DO NOT DISTURB ***** 

Execution can take from minutes to hours depending upon the size of the problem. The 

softkey menu returns to the default (pre-program) definitions. The remaining prompts occur 

after completion of the finite element calculations, and control the printing of output results. 

Computer prompt - Print Nodal Displacements ? (Y/N) 

User response - Print the nodal u, v, and w displacements, yes (Y) or no (N). Note the 

u displacement is not the total displacement, but only that part which is a 
function of y and z (i.e. U(y,z)). 

Computer prompt - Print Element Stresses and Strains ? (Y/N) 

User response - Print the element stresses, strains, strain energy density, and resultant 

nodal forces, yes (Y) or no (N). These quantities are computed and 
printed at each node of each element. 

The last line printed by FECAP for a normal execution is 

***** p ro biem Finished, Press RUN key for a New Problem ***** 


SAMPLE PROBLEM 

The problem selected for demonstrating the proper program input and output is a [0/90] 
symmetric graphite/epoxy plate exposed to a unifrom change in temperature. The region 
modeled by FECAP is shown in fig. 5. The proper user input is given below. User responses 
are followed by (ENTER) unless otherwise specified. 

DATA INPUT PHASE - Select a Key 
(Press softkey K0) 

Input Problem Title (<= 80 char.) 
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SAMPLE PROBLEM - [0/90] SYMMETRIC LAMINATE WITH THERMAL LOAD 
(ENTER) 

(Press softkey Kl) 

Input Number of Materials 

1 (ENTER) 

Is a Cylindrical Material Orientation Required ? (Y/N) 

N (ENTER) 

For Material No. 1 (Displayed at top of screen) 

Input Mat Prop's. Et,E2,E3,G23,G13,G12,NU23,NU13,NU12 
19.26E6,1.56E6,1.56E6,.487E6,.82E6,.82E6,.49,.238,.238 (ENTER) 

Input Thermal Prop’s. ALPHA 1 , ALPHA2, ALPHA3 
-.06E-6,15E-6,15E-6 (ENTER) 

Input Number of Angles 

2 (ENTER) 

Input Angles in degrees (one at a time) 

90 (ENTER) 

Input Angles in degrees (one at a time) 

0 (ENTER) 

DATA INPUT PHASE - Select a Key 
(Press softkey K2) 

Input Mesh Data File Name and MSUS (Ex. Fn,"700,0") 

MESH_4E,"700,0" 

Input the Maximum Allowable Semi_band width (Must be <= Array Dimensions) 

25 (ENTER) 

DATA INPUT PHASE - Select a Key 
(Press softkey K3) 

BOUNDARY CONDITION DATA - Select a Key 
(Press softkey Kl) 

NODAL BOUNDARY CONDITIONS - Select a Key 
(Press softkey K0) 

Specified Nodal Displacements: 

Input Direction of Specified Nodal Displacement(l-X,2-Y,3-Z) 

(Input 0 for End of Specified Nodal Displacements) 

1 (ENTER) 

Input Beg'ning node #, Ending node #, Increment, and Specified Value 
1,1, L0 (ENTER) 

Specified Nodal Displacements: 

Input Direction of Specified Nodal Displacement(l-X,2-Y,3-Z) 

(Input 0 for End of Specified Nodal Displacements) 

2 (ENTER) 

Input Beg'ning node #, Ending node #, Increment, and Specified Value 

1.7. 3.0 (ENTER) 

Specified Nodal Displacements: 

Input Direction of Specified Nodal Displacement! 1-X, 2- Y,3-Z) 

(Input 0 for End of Specified Nodal Displacements) 

3 (ENTER) 

Input Beg'ning node #, Ending node #, Increment, and Specified Value 

1.3. 1.0 (ENTER) 

Specified Nodal Displacements: 

Input Direction of Specified Nodal Displacement! 1-X, 2- Y,3-Z) 

(Input 0 for End of Specified Nodal Displacements) 

0 (ENTER) 

NODAL BOUNDARY CONDITIONS - Select a Key 
(Press softkey K2) 

Constrained Nodal Displacements: 
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Input Number of Different Constraint Conditions 
(Input 0 for No Constraint Conditions) 

1 (ENTER) 

Input Direction of Constrained Nodal Displacement(l-X,2-Y,3-Z) for Constraint 1 

2 (ENTER) 

Input Beg'ning node #, Ending node #, Increment, and Scale Factor 
3,9,3, 1 (ENTER) 

More Nodes for this Direction ? (Y/N) 

N (ENTER) 

More Nodes for this Constraint Condition ? (Y/N) 

N (ENTER) 

NODAL BOUNDARY CONDITIONS - Select a Key 
(Press softkey K4) 

DATA INPUT PHASE - Select a Key 
(Press softkey K4) 

GLOBAL LOAD CONDITIONS - Select a Key 
(Press softkey Kl) 

Input Average Axial Force and Temperature Change 
0,1 (ENTER) 

GLOBAL LOAD CONDITIONS - Select a Key 
(Press softkey K3) 

DATA INPUT PHASE - Select a Key 
(Press softkey K5 ) 

Input Printer Code for Output - CRT(l), PRINTER(701) 

701 (ENTER) 

Store Results on Output File ? (Y/N) 

N (ENTER) 

Print Nodal Coordinates ? (Y/N) 

Y (ENTER) 

Print Element Connectivity Data ? (Y/N) 

Y (ENTER) 

Print Boundary Condition Data ? (Y/N) 

Y (ENTER) 

Press Execute Softkey to Begin Execution 
(Press softkey K6) 

***** PROGRAM EXECUTING, DO NOT DISTURB ***** 

Print Nodal Displacements ? (Y/N) 

Y (ENTER) 

Print Element Stresses and Strains ? (Y/N) 

Y (ENTER) 

The output from the sample problem is given in Appendix C. It should be noted that the 
solution for this problem can be obtained from classical laminated plate theory (LPT) [5J. The 
stresses and strains predicted from FECAP are within 0.01 percent of the values predicted from 
LPT. 
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APPENDIX A 


Elements of fCbarl Matrix 

The specific form of the [Char] matrix is obtained using standard 4 ’th order tensor 
transformations [ 4 ], In contracted form [Cbar] can be written as a 6 x 6 symmetric matrix. The 
elements of [Cbar] for the plane orientation (fig. 3 ) are 

Cbarj i = m^Cj i + 2m 2 n 2 (Ci2 + 2C 66) + n 4 C22 
Cbar 12 = m 2 n 2 (Ci i + C22 - 4C66> + ( m4 + n4 ) c 12 
Cbar] 3 = m 2 Ci3 + n 2 C23 

Cbar 16 = mn(m 2 (Ci 1 - C12 - 2 C 66 ) + n 2 (Cj2 - C22 + 2 C 66 )) 

Cbar22 = n 4 Cj 1 + 2m 2 n 2 (Ci2 + 2 C 66 ) + m 4 C22 

Cbar23 = n 2 Ci3 + m 2 C23 (A.l) 

Cbar26 = mn(n 2 (Ci 1 - C12 - 2 C 66 ) + m 2 (Ci2 - C22 + 2 C 66 » 

Cbar 33 = C33 
Cbar 36 = mn(Ci3 - C23) 

Cbar44 = m 2 C44 + n 2 C55 
Cbar45 = mn(C55 - C44) 

Cbar55 = n 2 C44 + m 2 C55 

Cbar66 = m 2 n 2 (Ci 1 - 2C12 + C22) + (m 2 -n 2 )(m 2 -n 2 )C66 


where all other Cbarjj = 0, and m = cos(0) and n = sin(0). Note that Cbarjj = Cbarj, for 
i and j = 1 ,6. The elements of [C], the stiffness matrix in the principal material coordinate 
system, for an orthotropic material are given by [5] 

1 ■ v 23 v 32 v 12 + v 32 v 13 V13 + V12V23 

Cl 1 = c 12 = c 13 = 

E2E3Del EiE3Del EiE2Del 


1 ' v 1 3 V 3 1 v 23 + v 21 v 13 I-V12V21 

C22 = C23 = — - C33 = - — 

EjE3Del EiE2Del EiE2Del 

C44 = G23 C 55 = G 3 i Ci 6 = Gi2 

all other Cjj = 0 

where 


(A. 2) 


Del = 


1 - v 12 v 21 " v 23 v 32 - v 1 3 V 3 1 - 2v 21 v 32 v 13 
ElE 2 E 3 


(A.3) 



and Qj - Cji Ey, Gjj, and vy are the orthotropic Young’s moduli, shear moduli and 
Poisson's ratios, respectively, in the principal material coordinates. 

Note that for an orthotropic material there is a relationship between E and v which can be 
written as 


vy/Ej = Vjj/Ej (A. 4 ) 

The elements of [Cbar] for the cylindrical orientation (fig. 3 ) are written as 

Cbarj i = C'i ] 

Cbarj2 = m 2 C'i2 + n 2 C'i3 
Cbarj3 = n 2 C'i 2 + m 2 C'i3 

Cbari4 = mn(C'i2 * C’13) 

Cbar 15 = nC' 16 

Cbari6 = mC'i6 (A. 5 ) 

Cbar22 = m 4 C'22 + m 2 n 2 ( 2 C '23 + 4C44) + n^C'33 

Cbar23 = m 2 n 2 (C'22 - 4C44 + C33) + (m 4 +n 4 )C'23 

Cbar24 = m 3 n(C'22 - C'23 - 2C44) + mn 3 (C'23 + 2C44 - C'33) 

Cbar25 = m 2 n(C'26 ' 2 C 45 ) + n3c 36 

Cbar 26 = m 3 C '26 + mn 2 (2C'45 + c ' 36 ) 

Cbar33 = m 4 C'33 + m 2 n 2 ( 2 C '23 + 4C44) + n 4 C’22 

Cbar34 = m 3 n(C'23 + 2C44 - C'33) + mn 3 (C'22 ' C'23 * 2C44) 

Cbar35 = m 2 n(C '36 + 2C'45) + n 3 C'26 
Cbar 36 = m 3 C '36 + mn 2 (C*26 - 2 C' 45 ) 

Cbar 44 = m 4 C' 44 + m 2 n 2 (C' 2 2 - 2C’ 2 3 - 2 C’ 44 + C'33) + n 4 C’ 44 
Cbar 45 = m 3 C ’45 + nin 2 (C'26 * C'45 - C'36) 

Cbar 46 = m 2 n(C’26 - C'45 - C’36) + n 3 C'45 
Cbar55 = m 2 C'55 + n 2 C'66 
Cbar 56 = mn(C’66 - C'55) 

Cbar66 = m 2 C'66 + n 2 C'55 

where all other Chary = 0 , and m and n have the same definition as previously given. The 
elements of [C] are either the elements of [Cbar] given by equ. (A.l) when a plane orientation 
is used in combination with a cylindrical orientation, or the elements of [C] when the 
cylindrical orientation is used alone. For the later case certain elements of [C] used in equ. 
(A. 5 ) are zero. 

The elements of {a} are obtained using a transformation that is consistent with the definition of 
engineering shear strain, as apposed to tensor shear strain. The plane orientation results in an 
{a} with elements written as 



(A. 6) 


a x = m 2 cq + n 2 et 2 
a y = n 2 oq + m 2 ct 2 
a z = a 3 

a X y = 2mn(cq - 0 . 2 ) 

where all other a = 0, and m and n have the same definitions as previously given. 

The elements of {a} for the cylindrical orientation are written as 

«x = a 'x 

cty = m 2 a'y + n 2 a' z 
a z = n 2 a’ y + m 2 a' z (A.7) 

ay Z = 2mn(a'y - a' z ) 
a zx = na' X y 

a xy = ma 'xy 

where {a’} is equal to {a} from equ. (A.6) when a combination of plane and cylindrical 
orientations are used, and (a'} is defined as 

{a’ x , a’y, a' z , a’ yz , a' zx , a' xy ) = {otx ,<X2,ot 3 ,0,0,0} (A.9) 

when only a cylindrical orientation is used. 
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APPENDIX B 


Format for Mesh Data File 


The proper format for creating the mesh data file used in FECAP is given below. The 
statements are written in HP Basic 3.0. 

CREATE BDAT Fn_meshdat$&":,"&Msus_meshdat$,l, 
2*(Npe+2)*Nem+8*2*Nods+6 

ASSIGN @Filel TO Fn_meshdat$&":,"&Msus_meshdat$ 

OUTPUT @ FILE l;Nem, Nods, Npe 
FOR 1=1 TO Nods 
OUTPUT @ File 1;Y(I),Z(I) 

NEXT I 

FOR 1=1 TO Nem 

IF Npe=4 THEN OUTPUT @Filel;Nod(I,l),Nod(I,2),Nod(I,3), 
Nod(I,4),Imat(I),Iang(I) 

IF Npe=8 THEN OUTPUT @Filel,Nod(I,l),Nod(I,2),Nod(I,3), 
Nod(I,4),Nod(I,5),Nod(I,6),Nod(I,7),Nod(I,8),Imat(I),Iang(I) 

NEXT I 

ASSIGN @Filel TO* 


The variables in the above statements have the following definitions 


Fn_meshdat$ 

Msus_meshdat$ 

Nods 

Nem 

Npe 

Y(I),Z(I) 

Nod(I,J) 

Imat(I) 

lang(I) 


= File name of mesh data 
= Mass storage unit specifier 
= Number of nodes 
= Number of elements 
= Number of nodes per element, 4 or 8 
= Global nodal coordinates 

= Node numbers of I'th element, J=l,4 or 1,8 (Must be in consecutive 
CCW order around element) 

= Material number of I'th element 
= Angle number of material of I'th element 


Format for Boundary Condition Data File 

The boundary condition data file is created with the following format 
Mxnsl=MAX(Nsl(*)) 

CREATE BDAT Fn_bcdat$&":,"&Msus_bcdat$,l,8*(3’ ii Nsdf+3*Nsbf+ 
Ncd*3*Mxnsl+Ncd+3) 

ASSIGN @File3 TO Fn_bcdat$&":,"&Msus_bcdat$ 

OUTPUT @File3;Nsdf,Nsbf,Ncd 
FOR 1=1 TO Nsdf 

OUTPUT @File3;Ndd(I),Vbdf(I),Ibdf(I) 

NEXT I 

FOR 1=1 TO Nsbf 

OUTPUT @File3;Ndb(I),Vbsf(I),Ibsf(I) 

NEXT I 

FOR 1=1 TO Ned 
OUTPUT @File3;Nsl(I) 
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FOR J=1 TO Nsl(I) 

OUTPUT @File3;Nds(I,J),Ibcn(I,J),Vbcn(I,J) 
NEXT J 
NEXT I 

ASSIGN @File3 TO * 


The variables in the above statements have the following definitions 


Fn_bcdat$ 

Msus_bcdat$ 

Nsdf 

Nsbf 

Ned 

Ndd(I) 

Vbdf(I) 

Ibdf(I) 

Ndb(I) 

Vbsf(I) 

Ibsf(I) 

Nsl(I) 

Nds(I,J) 

Ibcn(I,J) 

Vbcn(I,J) 


= File name of boundary condition data 
= Mass storage unit specifier 
= Number of specified nodal displacements 
= Number of applied nodal forces 
= Number of nodal constraint conditions 
= Node numbers of specified nodal displacements 
= Magnitudes of specified nodal displacements 
= Directions of specified nodal displacements 
= Node numbers of applied nodal forces 
= Magnitudes of applied nodal forces 
= Directions of applied nodal forces 
= Number of nodes for each constraint condition 
= Node numbers for each constraint condition 
= Directions for each constraint condition 
= Scale factors for each constraint condition 


Format for Output Data File 

The output data file is created with the following format 

Max_len=INT(((Nods*3+Nem*(Npe*13+Npe*3)+3)*8+96)/256) 
CREATE BDAT Fn_outdat$&":,"&Msus_outdat$,Max_len+l,256 
ASSIGN @File2 TO Fn_outdat$&":,"&Msus_outdat$ 

OUTPUT @File2;T$,Nem,Nods,Npe 

FOR 1=1 TO Nods 

L=3*I-2 

M=3*I-l 

N=3*I 

OUTPUT @FiIe2;Gf(L),Gf(M),Gf(N) 

NEXT I 

FOR N=1 TONem 
FOR K=1 TO Npe 

OUTPUT @File2;sig(l),sig(2),sig(3),sig(4),sig(5),sig(6) 

OUTPUT @File2;Eps( 1 ),Eps(2),Eps(3),Eps(4),Eps(5),Eps(6), Energy 
NEXT K 

OUTPUT @File2;Fxyz(*) 

NEXT N 

OUTPUT @File2;Fxtot,Fytot,Fztot 
ASSIGN @FiIe2 to * 

The variables in the above statements have the following definitions 

Fn_outdat$ = File name of output data 

Msus_outdat$ = Mass storage unit specifier 

T$ = Problem title 

Gf(*) = Nodal displacements 

sig(*) = Element nodal stresses 
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Eps(*) = Element nodal strains 

Energy = Element nodal strain energy density 

Fxyz(*) = Element nodal resultant forces 

Fxtot = Summation of resultant x-direction nodal forces 

Fytot = Summation of resultant y-direction nodal forces 

Fztot = Summation of resultant z-direction nodal forces 

Those variables not defined above have the same definitions as given previously in this 
appendix. 
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APPENDIX C 


**«*4tt**44**»*44»*44*4*4«««4tf*«4*444*44««*«*>«***«*«44«4*44**4«4**44«***<***444 

Ft CAP PROOF AH 

27 Jul 1988 14:51:12 

444444444444444444444444444444444444444444444444444444444 4 4444444444444444 4 4444 

SAMPLE PROBLEM - [0/901 SYMMETRIC -AMI NATE WITH THERMAL LOAD 

44 4 4444444444444444444 4 444444444 4 44444 44444 44444444 444444^4444 4 4 44444 444 4 4 44 4 4 44 

INPUT DATA 

4 444 444444444444 444444 4444444444444444444444*4 4444444444 44444444444444 4- 444 4 44444 


CONTROL 1 N P U R MAT 1 G N : 


NUMBER OF ELEMENTS 4 

NUMBER OF NODES 9 

NUMBER OF NODES PER ELEMENT 4 

NUMBER OF IT A 1 E R I A L S 1 

MESH DATA FILE NAME MESH 4E 


BOUNDARY CONDITION DATA FILE NAME 
OUTPUT RESULTS FILE NAME 


MATERIAL PROPERTIES: 


MATERIAL NO. 1 

El 1.92 0 E + 0 7 

E2 1 . 5 6 0 E + 0 6 

E 3 . 1 . 5 6 0 E + 0 6 

G2i 4 . 870E+05 

613.. B.20QE+05 

012 8 . 200E+Q5 

NU 2 3 4 . 900E-0 1 

NU13 2. 380E-O1 

N U 1 2 2. 380E-01 

ALPHA 1 -6. OOOE-Oo 

AI.PHA2 1.5 OOF. -05 

ALPHA 3 1. 500E -05 


ANGLE INFORMATION FOR MAT. NO. 1 
ANGLE NO. THE TAT DEG.) 

1 9 0 . 0 0 

2 0.00 


NODAL COORDINATES: 

I'-CCOR 1)1 NATE 
0 . UOOE+OO 


NODE NO. 
1 
l 

4 

5 

6 


d . i.< 0 0 E - 0 1 
1 . OOOE+OO 
0 . OOOE+OO 
5. O00E-01 
1 . OOOE+OO 


/-COORDINATE 
0 . 0 0 0 E + 0 0 
0. OOOE+oO 
0 , OOOE+OO 
5. OOOE-'M 
5 . 0 0 0 E - o 1 
5 . OOOE-o 1 
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7 

8 
9 


ft . uftOE+ftft 
5. OOOE-Oi 
1 . OOOE+OO 


1 . OOOE+OO 
1 . OOOE+OO 
1 . OOOE+OO 


ELEMENT CONNECTIVITY DATA: 

ELEMENT NO. 

1 

2 

3 

4 


NODE NUMBERS 
12 5 4 

2 3 6 5 

4 5 8 7 

5 6 9 8 


MAT. NO. 
1 
1 
1 
i 


ANGLE NO. 
1 
1 
2 
2 


SPECIFIED DEGREES OF FREEDOM: 
NODE NUMBER DIRECTION 


1 

1 

4 

7 

1 

2 

3 


1 

2 

2 

2 

3 

3 


VALUE 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 
0. OOOE+OO 


CONSTRAINED NODAL DISPLACEMENTS: 

NODE NUMBER DIRECTION SCALE FACTOR 
3 2 1. OOOE+OO 

6 2 1. OOOE+OO 

9 2 1. OOOE+OO 

BANDWIDTH IS 21 

UNIFORM APPLIED AXIAL FORCE 0. OOOE+OO 

UNIFORM APPLIED DELTA TEMP 1. OOOE+OO 

UNIFORM COMPUTED AXIAL STRAIN 1.293E-06 


EXECUTION TIME IS 5.434E+01 SECONDS 

it*************:******************************#***************!******************** 

OUTPUT RESULTS 

**««t*************«*»«ft*4******ft««««««#**«*******«««»********«*«*«*««*********** 


NODAL DISPLACEMENTS: 


NODE NO. 
1 


4 

5 

6 

7 

8 
9 


U-DISPLACEMENT 
0. OOOOE+OO 
-1.51 65E-2 1 
-3. 0330E-21 
-7. 9 8 1 0 E — 2 2 
-1.51 65E-21 
-2. 2349E-21 
- 1 . 3369E-2 1 
-1.51 65E-21 
-1.6961E-21 


•DISPLACEMENT 
0. OOOOE+OO 
6. 4632E-07 
1 . 2926E-06 
0. OOOOE+OO 
6. 4632E-07 
1 . 2 9 2 6 E - 0 6 
0. OOOOE+OO 
6.4632E-07 
1 . 2926E-06 


W-UI SPLACEMENT 
0. OOOOE+OO 
0. OOOOE + OO 
0. OOOOE+OO 
1 . 0665E-05 
1 . 0665E-05 
1 . 0665E-05 
2. 1 329E-05 
2. 1 329E-05 
2. 1 329E-05 


ELEMENT STRESSES, STRAINS, STRAIN ENERGY, AND RESULTANT NODAL FORCES: 
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ELEMENT NO. 1. GLOBAL NODES 1 2 5 4 

Y COORD O.OOE+OO 5.00E-01 5.00E-01 O.OOE+OO 

Z COORD O.OOE+OO O.OOE+OO 5.00E-01 5.00E-01 

NODE EPSXX/SIGXX EPSYY/BIGYY EPSZZ/SIGZZ 

GAMMAYZ/TAUYZ GAMMAZX . TAUZX GAKMAXY/TAUXY 

1 1.293E-06/-2.098E+01 1.293E-06/ 2.09BE+01 2. 133E-05/ 2. 665E-14 

0. 000E+00/-3. 255E-32 - 1 . 596E-2 1 / -7 . 773E-16 -3.033E-21/ 8.194E-17 

SED = 1 . 580E-04 

2 1 - 293E-06/-2. 098E-+01 1.293E-06/ 2.098E+01 2. i 3 3 E - 0 5 / 5.329E-15 

0. 000E+00/-7 . 670E-48 -3. 762E-37/-1 . 832E-31 -3.033E-21/ 8. 194E-17 

SED = 1 . 580E-04 

5 1 . 293E-06/-2. 098E+0 i 1.293E-06/ 2.09BE+01 2.133E-05/ 5.329E-15 

-1.016E-20/-8.335E-15 -3. 762E-37/-3. 904E-31 -1 . 437E-21/ 1.391E-15 

SED = 1 . 5B0E-04 

4 1.293E-06/-2.09BE+01 1.293E-06/ 2.098E+01 2. 133E-05/ 2. 665E-14 

-1.016E-20/-8. 335E-15 -1 . 596E-21/-7. 773E-16 -1.437E-21/ 1.391E-15 

SED = 1 . 580E-04 

RESULTANT NODAL FORCES - 

NODE FX FY FZ 

1 4.129E-31 -5 . 244E+00 -5.059E-15 

2 1.943E-16 5. 244E+00 -4.711E-15 

5 1.739E-16 5. 244E+00 2.106E-1S 

4 -3 . 682E- 1 6 -5.244E+00 7.664E-15 

ELEMENT NO. 2. GLOBAL NODES 2 3 6 5 . 

Y COORD 5.00E-01 l.OOE+OO l.OOE+OO 5.00E-01 

Z COORD O.OOE+OO O.OOE+OO 5.00E-01 5.00E-01 

NODE EPSXX/SIGXX EPSYY/SIGYY EPSZZ/SIGZZ 

GAMMAYZ/TAUYZ GAMMAZ X , TAUZ X GAMMAXY/TAUXY 

2 1.293E-06/-2.098E+01 1.293E-06/ 2.098E+01 2. 1 3 3 E - 0 5 / 5. 329E-15 

0. 000E+00/-7 . 670E-48 -3. 762E-37/-1 . 832E-31 -3.033E-21/ 8.194E-17 

SED = 1 . 5B0E-04 

3 1 . 293E-06/-2. 098E+01 1.293E-06/ 2.098E+01 2.133E-05/ 1.954E-14 

0 . OOOE+OO/ 3 . 255E-32 1.596E-21/ 7.773E-16 -3. 033E-21 / B.194E-17 

SED = 1.580E-04 

6 1.293E-06/-2.098E+01 1.293E-06/ 2.098E+01 2.133E-05/ 1.954E-14 

6. 776E-21 / 5.557E-15 1 . S96E-21 / 7.773E-16 -1.437E-21/ 1.391E-15 

SED = 1 . 580E-04 

5 1 . 293E-06/-2. 098E + 01 1 . 2 9 3 E - 0 6 / 2.098E + 01 2.133E-05/ 5.329E-15 

6 . 7 76E-2 1 / 5.557E-15 -3. 762E-37/-4. 502E-32 - 1 . 437E-2 1 / 1.391E-15 

SED = 1 . 580E-04 

RESULTANT NODAL FORCES - 

NODE FX FY FZ 

2 -1.943E-16 -5.244E+00 -3 . 1 63E- 1 5 

3 1.726E-31 5.244E+00 -3.054E-15 

6 3.682E-16 5.244E+00 4.790E-15 

5 - 1 . 739E- i 6 -5.244E + 00 1.427E-15 


ELEMENT 

NO. 3, GLOBAL NODES 

4 5 

8 

7 


Y COORD 

O.OOE+OO 5.00E-0I 

5.00E-01 

O.OOE+OO 


Z COORD 

5.00E-01 5.00E-01 

1 . OOE+OO 

1 . OOE+OO 

NODE 

EPSXX/SIGXX 

EPSYY/SIGYY 

EPSZZ/SIGZZ 


GAMMAYZ/TAUYZ 

GAMMAZX. TAUZX 

GAMMAXY/TAUXY 


4 1 . 293E-06/ 2.098E+01 1 . 293E-06/-2. 09BE+01 2.133E-05/ 2.665E-14 


- 1 . 01 6E-20/-4 . 950E- 1 5 - 1 . 078E-2 1 / -3 . 837E- 1 6 -1 . 437E-21/-1 . 178E-15 

SED = 1 . 580E-04 

5 1.293E-06/ 2.09BE+01 1 . 293E-06/-2. 098E+01 2. 133E-05/ 5. 329E-15 

-1. 165E-20/-5. 672E-15 -1 . 505E-36/-1 . 234E-30 - 1 . 437E-2 1 /-I . 1 78E- 1 5 
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SED = 1 . 580E-04 

8 1 . 293E-06/ 2.098E+01 1 . 293E-06/-2. 098E+01 2.133E-05/ 5.329E-15 

-2.033E-20/-9.900E-15 - 1 . 505E-36/ - 1 . 234E-30 -3.592E-22/-2.946E-16 

SED = 1 . 580E-04 

7 1.293E-06/ 2.098E+01 1 . 293E-06/-2. 098E+01 2. 133E-05/ 2.665E-14 

-2.033E-20/-9. 900E-15 - 1 . 078E-21 /-8 . 837E- 1 6 -3. 592E-22/-2. 946E-16 

SED = 1 . 580E-04 

RESULTANT NODAL FORCES - 

NODE FX FY FZ 

4 3. 682E-1& 5.244E+00 -3.467E-15 

5 -1.473E-16 -5. 244E+00 -4.527E-15 

8 -2. 209E-16 -5 . 244E+00 9.237E-17 

7 -1.972E-31 5.244E+00 7.901E-15 

ELEMENT NO. 4. GLOBAL NODES 5 6 9 8 

Y COORD 5.00E-01 l.OOE+OO l.OOE+OO 5.00E-01 

Z COORD 5.00E-01 5.00E-01 l.OOE+OO l.OOE+OO 

NODE EPSXX/SIGXX EPSYY/SIGYY EPSZZ/SIGZZ 

GAMMAYZ/TAUYZ GAMMAZX , TAUZX GAMMAXY/TAUXY 

5 1.293E-06/ 2.098E + 01 1 . 293E-06/-2. 09BE+01 2. 133E-05/ 7. 105E-15 

5. 294E-2 1 / 2. 5 7 BE - 1 5 - 1 . 505E-36/-1 . 234E-30 - 1 . 437E-21 /-1 . 178E-1 5 

SED = 1.5B0E-04 

6 1.293E-06/ 2.09BE+01 1 . 293E-06/-2. 09BE+01 2. 133E-05/ 7. 105E-15 

8.470E-21/ 4.125E-15 1 . 078E-2 1 / 8 . B37E- 1 b -1.437E-21/-1. 17BE-15 

SED = 1 . 58QE-04 

9 1 . 293E-06 / 2. 098E + 01 1 . 293E-06/-2. 09BE + 01 2 . 1 3 3 E - 0 5 / 7. 105E-15 

6.776E-21/ 3. 300E-15 1.078E-21/ 8.837E-16 -3. 592E-22/-2. 946E-16 

SED = 1 . 580E-04 

8 1 . 293E-06/ 2. 098E+01 1 . 293E-06/-2. 098E+01 2. 133E-05/ 7.105E-15 

6 . 776E-2 1 / 3.300E-15 -1 . 505E-36/-1 . 234E-30 -3. 592E-22/-2. 946E-16 

SED = 1 . 580E-04 

RESULTANT NODAL FORCES - 
NODE FX FY FZ 

5 1.473E-16 5. 244E+00 -2. 801E-15 

b -3.6B2E-1S -5. 244E+00 -1.640E-15 

9 -2.650E-31 -5 . 244E+00 3.909E-15 

8 2. 2 0 9 E - 1 6 5 . 244E + 00 5.320E-16 


EQUILIBRIUM CHECK - SUM OF NODAL FORCES! 

Fx = 7 . 396E-32 Fy = 3.553E-15 Fz = -9.861E-32 

»«**»»»*****************»***•***»**«»•»«*•»*********«*****«*«****«******»******* 
»#*** Problem Finished, Press RUN Key Tor New Problem ***** 
******************************************************************************** 
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Table 1. Suggested System Requirements 


Computer 

o HP9000 Series 200 computer with HP Basic 3.0 operating system 
o 2Mbyte RAM 

o Floating point math co-processor card 
o Basic 3.0 compiler software 

Peripherals: 

o Dual 3.5 or 5.25 inch floppy disc drive 
o 15Mbyte hard disc 
o Dot-matrix printer 
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Figure 1. Coordinate system for generalized plane strain analysis of fiber reinforced 
composites. 
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Figure 3. Material orientation definitions. 
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Figure 4. Schematic diagram of stiffness matricies. 
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Figure 5. Sample problem geometry and FE mesh. 
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